HITTING DENSITIES FOR SPECTRALLY POSITIVE STABLE 

PROCESSES 



THOMAS SIMON 

Abstract. A multiplicative identity in law connecting the hitting times of completely 
asymmetric a— stable Levy processes in duality is established. In the spectrally positive 
case, this identity allows with an elementary argument to compute fractional moments 
and to get series representations for the density. We also prove that the hitting times are 
unimodal as soon as a < 3/2. Analogous results are obtained, in a much simplified manner, 
for the first passage time across a positive level. 



1. A MULTIPLICATIVE IDENTITY IN LAW 

Let {Xt, t > 0} be a spectrally positive Levy a— stable process (1 < a < 2), starting from 
zero and normalized such that 

(LI) E[e-^^*] = e*^", t, A > 0. 

For every x > 0, consider the first exit times 

= inf{t > 0, Xt > x}, = inf{t > 0, < -x} 

and the first hitting times 

Tx = mi{t > 0, Xf = x}, fx = inf{t > 0, Xt = —x}. 

Notice that since X has no negative jumps, one has = a.s. and is a positive stable 
random variable with index 1/a - see e.g. Theorem 46.3 in [H]. Besides, (II. ip and the 
optional sampling theorem give the normalization 

On the other hand, the process X crosses the level x > by a jump so that X^^ > x and 
Tx > Tx a.s. Introducing the overshoot = X^^ — x, the strong Markov property at time 
Tj. entails that {Xf = Xx^+t — ^t^, t > 0} is independent of {T^, K^). In particular, setting 
fy = inf {t > 0, Xf = —y} for y > 0, we have 

= + f^K. = + = x^T, + K'^h) 

where fi is an independent copy of fi and the two identities in law follow at once from the 
self-similarity relationships 

(Tj, T^, fx) = x"(Ti, Ti, fi) and = xKi, 
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which are themselves simple consequences of the self-similarity of X with index 1/a. Speci- 
fying the above to a; = 1 yields the basic identity 



(1.2) 



which has been known for a long time [6] . Notice that the laws of the three random variables 
appearing in the right-hand side of (11.2p are more or less explicit: the density of fi is that 
of a positive stable variable, the density of comes from a classical computation by 
Port involving some Beta integral - see e.g. Exercise VIII.3 in [2] - and the density of Ti 
can be obtained in changing the variable in the main result of [Ij thanks to the identity 
Ti = S'f", with the notation 5** = sup{Xs, s < t} for the supremum process. However, 
lack of explicit information on the law of the bivariate random variable {Ti,Ki) prevents 
from applying (11.21) directly to derive an expression of the law of ti. Recently Peskir [T2] 
circumvented this difficulty with an identity of the Chapman-Kolmogorov type linking the 
laws of Xi, Si and ri, and providing an explicit series representation at +oo for the density 
of Ti. In this paper we will follow the method of our previous article p[5j in order to derive, 
firstly, a simple identity in law for ti. Our main result reads: 

Theorem. One has 

(1.3) n = u^x fi, 

where U is an independent, positive random variable with density 

-(sin 7ra)t^/° 



fuAt) 

Proof. From (11. 2p we have for any A > 



7r(t^ — 2t cos na + 1) 



E 



E 

E [e 



E 



{Ti,K 



-{X°'Ti+XKi) 



] = E[, 



where in the third equality we recalled that fi is a (1/a)— stable positive variable with 
given normalization. On the other hand, the double Laplace transform of (Ti, Ki) had been 
evaluated long ago by Fristedt, and yields the simple expression 



a 



for every s > 1 - see (2.16) in [6]. It is possible to invert the right-hand side in noticing that 



1 



e-'^e-^dX 



and 



e"^^ (aA"-^E;(A")) dX 



for every s > 1, where E'^ is the derivative of the Mittag-Leffier function 



E 

n=Q 



r(l + an) 



z e 
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Above, the computation involving had been made originally by Humbert and we refer 
to the discussion after (5) in for more details. Putting everything together gives an 
expression of the Laplace transform of ri: 

(1.4) E[e-""^^] = Fi(A) - aF^(A), A > 0, 

where we used the notation Fa{x) = i?Q,(x") (notice that Fi{x) = e^) for every x > 0. We 
will now prove that the function Ga '■ A t— -Fi(A) — aF^(A) is completely monotone, reasoning 
as in Theorem 1 in |T5]. The well-known curvilinear representation of the analytic function 
1/r yields 

^i(^) = 77- / and «F^(A) = — / -——-dt 



for every A > 0, where H\ is a Hankel path encircling the disk centered at the origin with 
radius A - see e.g. Chapter IX. 4 in [1]. We deduce 



Ga{X) = ^ I e'gxit)dt, X > 



where 

/ 1 \ / aA"^i 



t-xj v^"-^", 

is analytic on C/{A} U (— oo,0]. Besides, since g\(t) —)•(« — 1)/2A as t — )■ A, it can be 
continued on the whole complex plane cut on the negative real axis, which entails 



/ e'gxit)dt, A > 



where if is a Hankel path independent of A. We can then compute, as in Theorem 1 in [15J, 
GaiX) = hm— / e~M - + - ^ - e"M r + — ^ ] ds 

-asinvra f°° / A°-^s" , , 

e — ^ ds 



TT /n V - 2S"A" COS Tia + A2° 



—a sm vra 



e-^' \ 1 ds 



Jo V ■s^" - 2s" COS 7r« + 1 

so that from (11. 4p and the expression of the density fu^ 



for all A > 0. Reasoning now exactly as in Theorem 3 in [15] entails finally the desired 
identity in law ti = Ua x ti. ^ 

Remark. In Theorem 1 of [6], the double Laplace transform of the function {t, x) i— )• F[tx > t] 
was computed for general Levy processes with no negative jumps. Our result can be viewed 
as an inversion of this Laplace transform in the particular stable case. 
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2. Fractional moments and series representations 

In this section we will obtain two expressions for the density of ri (a density which is 
readily known to exist from our main result, see also [TD] for the general context) in terms of 
convergent series, the first one being useful in the neighbourhood of and the second in the 
neighbourhood of +oo. Finding different asymptotic expansions for densities of first passage 
times is a classical issue in probability, see Chapter 1 in [8]. Recall also that this work has 
been done long ago for fi, respectively by Pollard and Linnik - see (14.31) and (14.35) in 
[Ti] . Our approach is different from the above references and relies on Mellin inversion. We 
will hence first compute the fractional moments of simple consequence of our main 

result: 

Corollary 1. For any s G (—1 — l/a;,l — l/a), one has 

^ ^ r(l - as) sin(7r(a - l)(g + l/a)) 
r(l - s) sin(7r(s + l/a)) 

Proof: From (11. 3p we have E [rf ] = E [ff ] E [U^] for every s G M. The computation of E [ff] 
is easy and classical: 

r(l - as) 



for every s < l/a. Similarly, one can show that 

-| poo 1 roo 

E[U^] = - -/ A-("^+i)G„(A)rfA = X-^-^^'\e'-aX--'E'^{X-))dX 

T{-as) Jo T{-as) Jq 

but I could not carry the computations further with the sole help of Mittag-Leffler functions. 
Instead, one can rest upon the residue theorem in order to evaluate directly 

-sinvra p x^+^Z" 

E f/^ = / dx. 

71 Jq x"' — 2a; cos 7ra + 1 

Setting A = s + l/a + 1 G (0, 2) and 



f{z) 



— 2z cos vra + 1 
one sees from Formula VIII. (10.2.3) in |1] that 

nK] = - (l^) (Res../ + ResU) 

where Xa = e'*-^""^'^ and Xa = e^^""^-''^ are the two conjugate roots of — 2x cos vra + 1. After 
some simple computations, one gets 

'sin(7r(A - l)(a - 1))A sin(7r(A - l)(a - 1)) sin(7r(a - l)(s + l/a)) 



E 

\^ sinvrA J sin7r(A — 1) sin(7r(s + l/a)) 

for every s G (—1 — l/a, 1 — l/a), which together with the previous computation for fi 
completes the proof (notice that 1 — l/a<l/a). ^ 
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Remark. The fractional moments of Ua and ti are finite over the same strip (— 1 — l/a, 1 — 
1/a), whereas those of fi are finite over the larger strip (— oo, 1/a). 

We will now derive series representations for the density /^^ of ti, inverting the fractional 
moments which were computed just before. Setting M{z) = ]E[rf ] for every z E C such that 
Re{z) G (— 1 — 1/a, 1 — l/a), one has 

M(iA) = E[e'^"^], A G M, 

where ai = logri is a real random variable with finite polynomial moments by Corollary 1, 
and which is absolutely continuous with continuous density f„^{x) = f^^{e^) over M. The 
Fourier inversion formula for /^^ entails then 

(2.1) /^^(a;) = _l_£M(it)x-'*rft, x>0, 

and we will start from this formula to obtain our series representation. The one at zero is 
particularly simple: 



Proposition 2. For every a; > one has 



^^l/a+n-l 



n> 



Proof: Suppose first x < 1. We compute the integral in (12.11) with the help of the usual 
contour Tji joining —R to R along the real axis and R to —R along a half-circle in the 
trigonometric orientation, centered at the origin. Because x < 1, the integral along the 
half-circle is easily seen to converge to when i? — )■ oo, so that it remains to consider the 
singularities of L{t) = M(it)x~'* inside the contour. By Corollary 1, the latter are located 
at i{n + 1/a), n > 0. After standard computations, one finds that 

/r(2 + an)sin(n7ra)x"+i/°\ / ax^+^Z" 

ReSi(„+i/„)L = -1 — ■ — = -1 



nT{l + l/a + n) J \T{-an)T{l/a + n) 

for any n > 0. By the residue theorem, this finishes the proof for x G (0, 1) (notice that 
ReSi/cfL = 0) and by analyticity, one can then extend the formula over the whole M"*". ^ 

The series representation at infinity is slightly more complicated than the one at zero, 
involving actually two series. Up to some painless normalization and after some simplifica- 
tions on the Gamma function, it had already been computed by Peskir [T2j with an entirely 
different, probabilistic argument involving a previous computation made for /^^ in [1]. 

Proposition 3. For every x > one has 

frAx) = J2Y(an)r(l/a-n) ^ ^r(-n/a)n\' 

Proof: By analyticity, it suffices to consider the case x > 1. We compute the integral in (12. ip 
with a contour Tn which is the reflection of Tji with respect to the real axis. Because x > 1, 
the integral along the half-circle converges to when R — )■ oo, so that one needs again to 
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consider the singularities of L inside r^^, which are located at i(l/Q; — n) and —in/ a, n >0. 
The term 

^r(an)r(l/«-n) 

follows then from the computations of Proposition 2 in replacing n by — n and taking into 
account the clockwise orientation of Tr. The last term follows from the simple computation 

R'eSi(„+i/Q,)L 



r{—n/a)n\ 

□ 

Remark. The two above representations are integrable term by term, which yields two series 
representations for the distribution function of ti : 

nn<x\ - 2^Y(-an)r(l + l/a + n) 

n>l \ / \ / / 

and 

. , „,™l/a-n ^ ^—n/a 

^^'-""^ = 1^ r(«n)r(l + l/a - n) " ^r(l-n/«)n!' 

One can also differentiate term by term in order to get expansions for the successive deriva- 
tives of : for any p > 1 one has 



. , ^„l/a+n-l-p 



„>-^r(-«^)r(i/a + n-p) 

and 



ri> 



T{an)T{ll a — n — v) T(—n/a — p)n\ 

1 ^ ' ^ ' ' n>l ^ ' ' 



The first term of all these expansions give the behaviour of the function at zero respectively 
at infinity. For instance, one has 



/ti(^) ~ T^/ „NT^/i , 1 /„N and f!^^{x) 



r(-a)r(l + l/a) ■'^1' ' T{-a)T{l/a) 

as x — >■ 0+, whereas 

^l/a-2 _Q,^l/a-3 

fri{x) ^ TTr I \ and f'{x) 



V{a-l)V{l/a) '^'^ ' V{(x)V{\l(x-2) 

as X — >■ +00. Notice finally that since a e (1,2), the function /^^ is ultimately completely 
monotone at infinity and /^^ ultimately completely monotone at zero. 
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3. Complements on the first exit time and the running supremum 

In this section, we will derive similar series representations for the density fj.^ of the first 
exit time Ti. As mentioned before, this random variable is connected to the supremum 
process {St, t > 0}, so that we will also get series representations for ^i. Our basic tool is an 
identity in law analogous to our main result, which was proved in our previous article |15j . 
Theorem 3. This identity reads 

(3.1) Ti = X fi 

where is an independent random variable with density given by 

^ _ -(sin7ra)(l + tV°) 
7ra;(t^ — 2t cos na + 1) 
over R+. As before, we begin in evaluating the fractional moments of Ti: 
Proposition 4. For every s G (—1, 1 — 1/a) one has 

^rysn ^ r(l + g)sin(7r/a) 
^ r(l + sa)sin(7r(s + l/a))' 

Proof. For every s G (— 1, 1 — 1/a) we have 

a 



1 /sin(7rs(« — 1)) sin(7rs(« — 1) — l/«) 



a \ sinvrs sin(7r(s + 1/a)) 

sin(7rsa;) sin(7r/a) 

a sin(7rs) sin(7r(s + l/a)) 

where the last equality follows from tedious trigonometric transformations. From (13. ip . the 
previous computation made for E [f*], and the complement formula for the Gamma function, 
we get 

sin(7r/Q;) /r(l — sa) sin(7rsa)) \ r(l + s) sin(7r/a) 



□ 



^ sin(7r(s + 1/a) V «r(l - s)sin(7rs) J 1(1 + sa) sin(7r(s + 1/a)) 

as desired. 



Remark. Similarly as above, the fractional moments of and ti are finite over the same 
strip (—1, 1 — I /a), whereas those of fi are finite over the larger strip (— oo, 1/a). 

We notice in passing that this computation allows to show that the law of the independent 
quotient Ti/Ti is Pareto(l — 1/a), a fact which had been proved in [S] in the general context 
of stable processes in duality: 

Corollary 5. The density of the independent quotient Ti/Ti is given by 

sm{n/a) 
7rii-i/"(l + t) 

over M"*". 
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Proof. By Mellin inversion, it is enough to identify the fractional moments. Proposition 4 
and the classical computation for the moments of Ti entail 



E 



fs 



sin(7r/a;) sin(7r/a)t'* 



sin(7r(s + 1/a)) Jq 7rti-V"(l + t) 



for every s G (— 1 — 1/a), where the second equality follows from Formula VIII. (10.2.4) 
in H. 

We now come back to series representations for frp_^, which are obtained from Proposition 
4 and Mellin inversion similarly as in Propositions 2 and 3. Since we choose exactly the 
same contours, we will leave all the details to the reader. Notice that contrary to ri, the 
two-series representation for Ti is the one at zero and that the first series therein defines an 
entire function. 

Corollary 6. One has both representations 

frM) = y2~Ff ^\■r^^ , i / V' 

^ ^-^ aT(an — IjFil + l a — n) 

n>l ^ / \ I I 

and 

n>l ^ ' n>l ^ / \ / / 

Thanks to the aforementioned identity Ti = 5*]"°, changing the variable gives correspond- 
ing series representations for the density /^^ = ax'^"'^^^ fj,^{x~°'). The first one had been 
obtained originally in [T] after solving some fractional integral equation of the Abel type. 
See also [TT] for a new proof involving Wright's hypergeometric function. We stress that up 
to the Wiener- Hopf factorisation which is a fundamental tool for this type of questions, over- 
all our argument to get this series only makes use of undergraduate mathematics. Besides, 
this simple method allows to get the series representation at infinity. 



Corollary 7. One has both representations 



fsA-) - gr(an-l)r(l + l/«-n)' 
and 

fa (x) = > z + > — : ; ; -, X > 0. 

n>l ^ ' n>\ ^ ) \ I ' ) 

Remark. Notice that one really needs the two series to get the behaviour of /j.^ and its 
derivatives at zero. For instance one has 

1 1 

/Ti(a;) r and ~ — , , , . /^^(x) 
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as X — 7- 0+. On the other hand, only one series is important to get the behaviour of f^^ and 
its derivatives at infinity: one has 



^-(a+2+p) 

f^iAx) ~ —, r and A (x) ~ — — - 



as X — >■ +00 for every p > 1 (this is in accordance with the general asymptotics derived in 
[7], see the final remark therein). 

4. Some remarks on unimodality 

Recall that a real random variable X is said to be unimodal if there exists a G M such 
that the functions F[X < x] and F[X > x] are convex respectively in (—00, a) and (a, +00). 
If X has a density fx, this means that fx increases on (—00, a] and decreases on [a, +00). 
We refer e.g. to Section 52 in [14] for more on this topic. Differentiating the density of Ua '■ 

we find that there is a unique positive root, so that the random variable Ua is unimodal. Since 
the same property is known to hold for fi - see e.g. Theorem 53.1 in [H] for a much stronger 
result, in view of our main result it is natural to ask whether the variable ti is unimodal as 
well. This seems also plausible from the behaviour of /r^ and /^^ at zero. In general, the 
product of two positive independent unimodal random variables is not necessarily unimodal. 
However the notion of multiplicative strong unimodality which had been introduced in [3], 
will allow us to give a quick positive answer in half of the cases: 

Proposition 8. For every a < 3/2, the random variable Ti is unimodal. 

Proof. Recalling that fi is unimodal for every a G (1,2), from (11. 3p and Theorem 3.6 in [3j 
is is enough to show that the function 1 1— )■ fuai^^) is log-concave over M, in other words that 
the function ga{t) = log(e* — 2cos7ra + e~*) is convex over M. Differentiating twice yields 

„ 4(1 — cos vra cosh t) 

[e^ — 2 cosTTo; + e ^j^ 

and we see that Qa is convex over M as soon as a < 3/2. ^ 

When a G (3/2, 2) the above proof shows that Ua is not multiplicative strongly unimodal 
(MSU) in the terminology of [3]. In order to get the unimodality of ti when a G (3/2,2), 
one could of course be tempted to obtain the MSU property for positive stable laws with 
index (3 > 2/3. Unfortunately, we showed in ^16j that this property does not hold true as 
soon as /3 > 1/2. Nevertheless, in view of the unimodality of fi and that of ri for a = 2 - see 
also Rosler's result for general diffusions on the line [13], we conjecture that ti is unimodal 
for every value of a. At first sight the problem does not seem easy because of the unstability 
of the unimodality property under multiplicative convolution. See however our final remark. 
Differentiating now the density of Ta : 

fT (t) = —r^ -,r((l-2a)t2-2at2"^/" + 2t(a-l)cos7ra + 2at^-^"cos7ra + l), 

Q;(t"' — 2tcos7ra; + 1)^ 
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we find after some simple analysis tliat tfiere is also a unique positive root, so tliat tlie 
random variable is unimodal. But contrary to Ua, one can clieck tliat it is never MSU, 
so that the above simple argument cannot be applied. Nevertheless we have the 

Proposition 9. For every a < 3/2, the random variable Ti is unimodal. 

Proof. From (13.11) and the so-called Kanter representation for positive stable laws - see 
Corollary 4.1 in [9j, we have the identity 

Ti = T„ X L^-" X b^{U) 

where L is a standard exponential variable, U an independent uniform variable over [0,7r], 
and ba{u) = (sin((l — 1 / a)u) / smi^u))"'^ sm{u / a) / sm{u) is a strictly increasing function 
from (0, 7i) to M^. Besides, from the beginning of the proof of Theorem 4.1 in [9], we know 
that the logarithmic derivative of ba is positive and strictly increasing, so that the same 
holds for b'^ itself, because ba is also positive and strictly increasing. By Lemma 4.2 in [9] 
we conclude that baiU) is unimodal: again, from (13. ip and Theorem 3.6 in [3], it suffices to 
show that the density t ^ fxcxL^-'^i^^) log-concave over M. We compute 



— (sin7ra)xi-« f°° i mq-i(1 + u , 

It^xl^-^x) = 7 7T— / exp-{x u)^- — — —du, x>0. 

7Ta{a — 1) Jq — 2u cos na + 1 

It is easily seen that u ^-^ l + u^^°' is log-concave over and (t, u) i— e*/u log-convex, hence 
convex, over M x M+. In particular, the function 



1 

{t,u) ^ exp - (eV^) (1 + 



u 



is log-concave over R x M"*" and by Prekopa's theorem, it is enough to show that the function 
u H-> ■u«^/('u^ — 2-ucos7ra + l) is log-concave over M"''. We compute its logarithmic derivative, 
which equals 

(2a — 3)u* + 4m^(2 — a) costto — (4(2 — a) cos^vra + 2a)u^ + AucosTia — 1 

{a — 1)^2(^2 — 2u cos 7Ta + 1)^ 

and we see that it is negative over R"*" as soon as a < 3/2. ^ 

Final Remark. The above proof does not work to obtain the unimodality of Ti when 
a > 3/2, and we do not know as yet how to tackle this situation. However, we can use the 
same argument to obtain the unimodality of ti when 3/2 < a < 1 + l/\/2 < 2. Details, 
which are technical, can be obtained upon request. All these questions will be the matter of 
future research. 
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